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Abstract: 

A micromachined deformable mirror with pixelated electrostatic 
actuators is proposed. The paper begins with a physical description of 
the proposed mirror. Then a mathematical model in the form of a nonlinear 
partial differential equation describing the mirror surface deformations 
is derived. This model is used to derive the required voltages for the 
actuators to achieve a specified static deformation of the mirror 
surface * This is followed by the derivation of a static nonlinear 
feedback controller for achieving noninteracting actuation. Then the 
structure for a complete control system for wavefront correction is 
proposed. The paper concludes with a discussion of the physical 
implementation of the proposed control system. 

1. INTRODUCTION 

In the development of large space interferometers and multi-aperture 
reflectors, deformable mirrors are used to compensate for distortions in 
elements of the optical train and/or in the instrument’s field of view. 
Such mirrors should be small and lightweight. Moreover, they should be 
highly pixelated so that the deformations can be controlled with high 
lateral resolution. 

The first actively controlled deformable mirrors were developed by 
NASA in the I960’ s for use as solar collectors or as ground-based 
telescopes [1]. Since then, there has been extensive development in this 
area. Comprehensive surveys of works on actively controlled deformable 
mirrors were given by Ealey [2] and Tyson [3]. In 1977, Grosso and Yell in 
[4] developed a membrane mirror whose deformations are controlled by 
means of discrete electrostatic actuators. Subsequently, various forms of 
deformable mirrors with discrete piezoelectric and magnetostrictive 
actuators were also developed [5] -[9]. The advent of silicon VLSI 
technology has made possible the integration of deformable mirrors with 
microelectronic circuitry. In 1983, Hornbeck [10] perfected a deformable 
mirror device with pixelated mirror elements whose size is 51 pm square. 
The mirror deformations are controlled by electrostatic actuators driven 
by microelectronic circuits which are integrated with the mirror 
assembly. His mirror was used primarily as a light modulator. Later, in 
1989, a wavefront control device with a deformable mirror integrated with 
control and sensor units was introduced by Ealey and Wheeler [11]. In 
their device, the actuators are spaced 1.0 mm apart. The voltages applied 
to the actuators are on the order of 200 volts. Recently, efforts have 
been initiated at the Jet Propulsion Laboratory in exploiting micro- 
machining technology to develop deformable mirrors with the afore 
mentioned characteristics and with pixelated electrostatic actuators 
which are spaced less than 25 microns apart. In this paper, attention is 
focused on the analytical design of control systems for such deformable 
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mirrors. 

The paper begins with a physical description of the proposed mirror. 
Then a mathematical model in the form of a nonlinear partial differential 
equation describing the mirror surface deformations is derived. This 
equation is used to derive the required actuator voltages to achieve a 
specified static deformation of the mirror surface. This is followed by 
the derivation of a static nonlinear feedback controller for achieving 
noninteracting actuation. Then the structure for a complete control 
system for wavefront correction is proposed. The paper concludes with a 
discussion of the physical implementation of the proposed control system. 

2. PHYSICAL DESCRIPTION OF DEFORMABLE MIRROR 

Figures 1 and 2 show respectively the sketches of the top and side 
views of the proposed deformable mirror with pixelated capacitive 
actuators. The mirror may be realized as "flip chip"-type assemblies 
consisting of two matched micromachined silicon structures mounted 
face-to-face and fused together along their peripheries. The key elements 
of the mirror consist of simple, easily replicated, electrostatic linear 
actuators, each responsible for pulling on a small portion of a thin 
flexible silicon membrane which is the substrate for the deformable 
mirror. The mirror surface is formed by depositing a metallic or 
multi-layer dielectric film on the membrane. The membrane with posts (See 
Fig. 2) is micromachined from a silicon sheet. The posts serve as supports 
for the membrane and also as halves of the electrostatic actuators. The 
bottom half of the mirror assembly consists of a set of posts with four 
silicon blades attached to each post. These blades serve as leaf springs 
for supporting the posts of the upper mirror assembly, and for providing 
a restoring force for the actuation system. This bottom assembly Is 
micromachined from a silicon wafer. The electrostatic actuators are 
formed by attaching conductive pads to the upper posts and the bottom 
half of the mirror assembly. The electronic element access, electronic 
actuator drivers, and possibly the feedback controller circuitry may be 
monolithically integrated into the mirror assembly. 

We note that the geometric structure of the deformable mirror proposed 
here differs from that of Hornbeck [10]. In his mirror, each actuator, 
when activated, produces a concave deformation of the mirror surface over 
the entire pixel. This causes focusing of the incoming light beam in 
front of the pixel. Here, each actuator pulls down on the mirror surface 
at a post area and thereby induces deformation over adjacent portions of 
the mirror surface. Except for the flat spots over the post areas, the 
overall shape of the mirror surface is determined by the displacements of 
all the actuators. 

The initial performance goals for the proposed mirror will be the 
control of a 32 x 32 pixel flat mirror with 10 nm accuracy. Once these 
goals have been achieved, efforts will be directed at extending the 
number of pixels/control elements until 10 nm accuracy can be achieved 
over a 1024 x 1024 pixel surface. 

3. MATHEMATICAL MODEL 

2 

Let Q be an open connected subset of the Euclidean plane R with a 
piecewise smooth boundary 5Q representing the spatial domain of the 
mirror. We introduce a mesh on Q whose mesh points are denoted by x = 

mn 

(x ,x ), m = 1 , . . . , M; n = For a rectangular mirror, Q is 

lm 2n 
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specified by £2 r = {(x^x^ e R 2 : |xj s I x g I * M, where the s 

are specified lengths. For a circular mirror, £1 is specified by the disk 
£1 = {( r 0), Osrsr.Oses 2n}. At each mesh point x , we 

introduce a patch £3 , a bounded open subset of £3 representing the 

mn 

effective spatial domain of the (m,n)-th actuator force containing as 

an interior point. Typical meshes and patches for the rectangular and 
circular mirrors are shown in Fig. 3. 

Let the mirror surface be a thin membrane with density p = p(x) being 
a specified positive piecewise smooth function satisfying the following 
bounds: 

0 < p — p ( x ) s p < +co for all x € £2. (1) 

min max 

The variation of the mass density due to the supporting posts can be 
included in p by setting p(x) = p (a known constant) for x € £2 ■ 

mn 

Let <r , i,j = 1,2 denote the components of the symmetric stress 

tensor in the mirror surface satisfying the positivity condition 


2 2 

C ll^ll 2 s y V (r £ £ SC ll£ll 2 for all £ = (£ J e r2 > ^ 

1 ifrijfri iJ 1 J 2 

where c and c are known positive constants. In the special case with 
1 2 

uniform tension T, we have <r = T5^, where 5^ denotes the Kronecker 
delta. 

The downward displacement u(t,x) normal to the mirror surface at a 
point x e £2 and time t £ 0 can be described by the following equation: 


p(x) 




(3) 


where f = f(t,x) is the surface force density whose explicit form will be 
derived later. Assuming that the mirror is rigidly attached to its 
boundary 3£2, u must satisfy the boundary condition 


u(t,x) = 0 for x e 3£2 and t £ 0. 


(4) 


Finally, the initial conditions for u are specified by 

u(0,x) = u (x), -lr(t.x) = u* (x) for x e £2. (5) 

o O t o 

To derive an explicit expression for the surface force density f, we 
first consider the electrostatic force density over a patch £2^ due to a 

specified voltage V (t) applied to the (m.n)-th actuator. We assume that 

the mirror surface curvature is small so that each patch £2^ is 

essentially parallel to the bottom assembly. Thus each actuator’s 
conductive surfaces can be regarded as making up a parallel-plate 
capacitor. Neglecting fringing effects of the electric field at the 
boundary of £2 , the electrostatic force density is given by 

mn 
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for all x e £1 


( 6 ) 


V (t) 


i ( * v > \ 

f f t > *■ mn 

(t,X) = C -= 7”7 7 - 

e 2 o D-u(t,x) 

' mn 


where D is the distance between the undeformed mirror surface and the 
bottom plane and e is the permittivity of free space. When D » 

o 

|u(t,x )|, (6) can be approximated by 


f (t,x) = c V 2 (t)/D 2 for all x € £2 . 

c Z o mn nn 


( 6 ’ ) 


Considering each leaf spring as a small cantilever beam having uniform 
cross section with moment of inertia I and Young’s modulus E, the force 
density f due to four leaf springs over the patch £2 is given by 


f (t,x) = 


12EI 
l 3 A 


u(t,x ) for all x e £2 , 


(7) 


where A denotes the area of the patch £2 . Here, we have neglected the 

mn mn 

inertial effects of the leaf springs. 

Let <p denote the spatial weighting function associated with the 

mn 

(m,n)-th actuator such that <p (x) = 0 for x e Q - {Q v dQ ) . 

mn mn m n 

Combining (6) and (7), equation (3) becomes a nonlinear partial 
differential equation given by 

2 2 


p(x) 


d 2 u 


at 


-[ E IT ('./■> &■) 

i = 1 j = l i v y 


M N 


•E E{t*.( 


V (t) 

mn 




D - u(t,x ) 

mn ' 


n=ln=l 


12EI 

i 3 a 


} 


• u(t, x ) V <f> (x) . 


( 8 ) 


Let K(x,x* ,t,x) denote the Green's function corresponding to the 
solution of the linear equation: 


om-tt -l i icr h ) ,x) Jr ) * a(t 

at i = i j=i i v y 


-t,x - x' ), 


(9) 


with boundary and initial conditions given by (4) and (5), where 5 
denotes the Dirac delta function at t = t and x = x’ . Equation (8) can be 
reformulated as a nonlinear integral equation: 

u(t,x) = f K(x,x’,t,0)u (x* ) dx f + f x ’ * t* 0)u’ ) dx’ 

i Q ° in dt ° 


Q 


M N „ 

r t r f 1 , V (t) 

> M n K<X ' X ’' t ' T) l I{t £ . D-“tx ) - 


m=ln=l 


12EI 
l 3 A 


u(t,x ) 

mn 


} <t> (X’ 

I mn 


)dx’ 


(10) 


Under the assumption that_ the mirror deformations ^ver £2 are 

mn 

sufficiently small compared to D so that (6) may be approximated by (6*), 
equations (8) and (10) become linear. This assumption may not be 
justified when D is made small so as to reduce the operating voltage 
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levels of the actuators. For example, in Hornbeck’ s deformable mirror, D 
is 620 nm, and the peak mirror deformation for normal operation is around 
100 nm. Evidently, (6’ ) is not a good approximation for this case. 


4. STATIC SHAPE CONTROL 


Let u = u (x) be the desired static shape of the mirror surface 
d d 

defined over the entire spatial domain 0. It Is required to determine the 

static voltages V for each actuator to achieve the desired shape u . 

mn a 

Let K = K (x, x’ ) denote the Green’s function associated with the 


boundary-value problem: 



with boundary condition 


(x) 



5(x - x’ ) , X € n, 


( 11 ) 


u(x) = 0 for x € dQ. 


( 12 ) 


Then, the static equation corresponding to (8) can be reformulated as 
an integral equation for the mirror displacement u corresponding to given 

static actuator voltages V : 

mn 


u(x) 



(x, x’ ) 


i [{-H 


m = 1 n = 1 


V 

mn 

D - u(x ) 

mn 



12EI 

A 


u(x 

mn 


’} 


mn mn 


<t> (x’ )dx’ , 

mn 


for X € Q. (13) 

Let u = (u(x ti ) u(x 1N } u(x mi ), . . . ,u(x mh ) )T and V k = (V^ 

V k V k ,...,V k ) T , k = 1,2. Setting x = x in (13) leads to a set of 

IN Ml MN ij 2 

M x N algebraic equations relating u and V . These equations can be 
written as 

(I + K)u = P(u)V 2 , (14) 


where I denotes the MN x MN identity matrix; K is the MN x MN matrix 


whose k-th row K is given by 
k 


= (K^ t . . . . . . ,K*J, . . . ,K^), i = integer (k/M), j = k Mod(M) , (15) 


where 


= i 
mn Uq 


K J = \ K (x , x' ) <p (x* ) dx 


s i j 


•} 


12EI 
l 3 A 


(16) 


P(u) is the MN x MN matrix whose k-th row is (p 1 ^ (u) , . . . , p 1 ^ (u) , . . . , 

11 In 


p^(u), p-( u )), Where 


Mi 


p 1J (u) = 

mn 


2 (D 


-jf 

- u(x ) ) 2 


K (x , x’ ) 4> (x* ) dx 

s 1 j mn 




i = integer(k/M) , j = k Mod(M). (17) 
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Explicit expressions for K and p lJ corresponding to rectangular and 

r s inn 

circular mirrors are given in the Appendix. 

If we set u = u d = u d (x i h } ‘ ' ’ U d (x mi } ’ ' ' ’ U d (X MN }) ’ 

then (14) becomes a set of M x N linear algebraic equations for the 

— “2 
unknown actuator voltages V. Evidently, if P(u^) is nonsingular, then V 

is uniquely determined by 

v 2 = P(u ) _1 ( I + K)u (18) 

d d 

The matrix P(u ) is nonsingular if and only if its rows are linearly 

d 

independent, or equivalently the M x N matrices given by 


P (u ) = 

i j d 


r P iJ (u ) • • 


*11 d 

• 

IN d 

• 

* 

• 

p 1 J (u ) • • 

# 

• p‘i(u ) 

*M1 d 

MN d 

independent. 

Since p 1 ^ 


i = 1.....M; j = 1 N, 


(19) 


displacement u , the set X> of all displacements u such that P(u d ) is 

nonsingular corresponds to the set of all mirror surface displacements 
which have one-to-one correspondence with the actuator voltages In 

fact, if we define the nonlinear mapping u Af(u) by N(u)=P(u) 1 ( I + K)u, 

then N is an invertible mapping with domain D. In the case where the 

desired mirror surface deformation u = u (x) has spatial symmetry, the 

d d 


number of equations in (14) can be reduced accordingly. 

Now, given a set of actuator voltages V, the corresponding mirror 
surface displacements at the mesh points x ^ can be determined by solving 


the nonlinear equation (14) for u. If we define the nonlinear mapping 
u A r (V 2 )u by N {y 2 )u = P(u)V 2 - Ku, then the solutions correspond to the 


fixed points of A'(V). In physical situations where the actuator voltages 
satisfy a magnitude constraint of the form V - V < the set of all 
admissible V’ s is given by V = (V e R MN :V 2 n ^ ^ ax ’ m = n = 

1 N} (a hypercube in R MN ). Then the set of all admissible u^s is 

given by N 1 (V ) . 


Remarks : 

(r-1 ) Once the required actuator voltages V for achieving the 

desired static mirror surface displacements at the mesh points are 
determined, the mirror surface displacements at other points in the 
spatial domain Q can be found using (13) with u ^ x mn ^ se ^ U d^ X mn^* 

(r-2) In all the existing works on deformable mirrors, it is 
assumed that the static mirror surface displacements at the mesh points 
x are related to the actuator inputs by a linear transformation 

mn 

commonly called the influence function matrix which is valid for small 
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displacements. For large mirror_ displacements, the relation between 

u(x ) and the actuator voltages V is given by (13) with x set to x . 
i j mn 1 J 

This relation is implicit and nonlinear. In the special case of small 

displacements such that approximation (6’ ) holds, this relation becomes 

linear. Using (13) with a single actuator at 12 to obtain u(x) in terms 

mn 

of V gives the usual influence function. 

mn 

(R-3) If the number of actuators is less than the number of mesh 
points at which the desired mirror displacements are specified, then, 
under the linear approximation (6’ ), (14) consists of a set _ 2 of 

overdetermined linear algebraic equations for the unknown variables V . 

We may use the least-squares solution, which corresponds to obtaining the 
pseudo-inverse of the influence function matrix [3]. 

Noninteracting Actuation: 

Due to stress in the membrane, the voltage V applied to the (m,n)-th 

mn 

actuator will influence the membrane displacements at all the patch 
locations. To simplify the mirror deformation control, it is desirable to 
introduce appropriate feedback and new control variables c such that 

c only influences the mirror displacement u at the mesh point x , and 

mn mn 

u (x ) only depends on c . Thus, the new controls c produce 

mn mn mn 

noninteracting actuation of the mirror surface displacements at the mesh 
points. To achieve noninteraction, we introduce a static feedback control 
of the form V 2 = F(c - u), where F is a feedback gain matrix which may 
depend on u. Thus, in view of (14), we have 

(I + K)u = P(u)F(c - u), (20) 

or 

[I + (I + K) -1 P(u)F]u = (I + K) -1 P(u)Fc, (21) 

where c = (c c c c ) T . To achieve noninteracting 

11 IN Ml MN j 

actuation, we seek an MN x MN matrix F such that (I + K) P(u)F = A, a 

specified constant diagonal matrix with nonzero diagonal elements A^. 

Thus, u and c are related by a diagonal matrix operator given by 

u = [I + A] ^c, (22) 

and the required feedback gain matrix F is given by 

F = P(u) _1 (I + K)A. (23) 

In physical terms, the static feedback control in effect cancels the 
spring coupling forces to produce noninteracting actuation with respect 
to the new control c. In the special case where the mirror displacements 
are small as compared to the actuator gap D so that approximation (6’ ) is 
valid, P becomes a constant matrix. Consequently, F is also a constant 
matrix and the noninteracting controller is linear. We note that P(u) and 
K depend on the parameters E, I, D, and t whose values can be 

accurately estimated. Therefore F can be determined with good accuracy. 
Finally, the foregoing noninteracting controller is also valid for the 
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dynamic case, since the couplings between the mirror displacements at the 
mesh points involve elastic forces only. 


5. DYNAMIC SHAPE CONTROL 

Assuming that the desired nominal static-shape for the mirror is 

attainable by appropriate choices of the actuator voltages V ' , it is of 

interest to control the deviations of the mirror shape about u d for 

wavefront correction. We propose to achieve this objective in three 
steps. Figure 4 shows the proposed structure of the overall control 
system for wavefront correction. First, a static feedback controller for 
achieving the desired nominal static mirror-shape and noninteracting 
actuation of the deformable mirror is introduced. The function of the 
minor-loop feedback controller is to modify the dynamic response of the 
actuators. This modification can also be performed ahead of the static 
feedback controller depending on the method of implementation. Finally, a 
global feedback controller which makes use of the output of the wavefront 
estimator to generate the appropriate actuating signals for wavefront 
correction Is Introduced. The static feedback controller has already been 
discussed in Sec. 4. In what follows, the discussion will be devoted to 
the modification of the dynamic response of the actuator, and the global 
shape controller for wavefront correction. 

5.1 Actuator Dynamics Modification 

The main objective here is to modify the dynamics of the actuator to 
ensure satisfactory response to input commands. Since this controller is 
to be integrated with the mirror assembly, the control law should have 
the following properties: 

(i) It should be sufficiently simple so that it can be realized by 
microelectronic circuitry which can be integrated monolithical ly with the 
mirror microstructure. 

(ii) It should be model independent so that it is unnecessary to 
identify the system parameters for their implementation. 

(iii) Its performance should be sufficiently robust with respect to 
system parameter variations. 

Since the mirror surface has very little internal damping, the actuator 
forces may induce undesirable surface vibrations. Therefore, it is 
necessary to introduce damping externally. A possible approach is to 
introduce external passive damping. This may be achieved by housing the 
bottom mirror assembly in an enclosure which contains air, and has minute 
holes for air passage. Alternatively, damping can be achieved by means of 
active feedback controls. 


To derive appropriate forms for the active feedback controls, we make 
use of the partial differential equation (8) linearized about the nominal 

A A 

mirror surface deformation u = u (x) produced by the static actuator 

d d ^ 

voltage V which produces the desired u^. Note that u^(x) = u d ^ x ^ on ^V a t 

A _ 

the mesh points x = x Let Su = u - u , and 5V = V - V. The linearized 

mn d 

equation (8) is given by 


p(x) 


3 2 5u 


at 


-l J -Sr ('./*> #) 

i = 1 j = 1 l K V 
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(24) 


c V 2 19T7T , e V 5V (t) 

12EI_] Su( t,x ) +• ° 

(D - u (x )) 3 l 3 A J mn (D 

ra = l n = 1 d mn mn d mn 


n n 

•H{(- 


V SV ( t ) \ 

:> mn mn I 

- u (x )) 2 f 


<P (x), 

mn 


with boundary condition 

<5u(t,x) = 0 for x e and t i 0, 
and initial conditions: 


3<5u 


5u(0 f x) = u (x) - u (x), a . 

o d Ot 


(0,x) = u’ (x) for x e n. 


(25) 


(26) 


about u^ given by 


Consider the total energy functional of the perturbed mirror surface 

-TtHSf '[(«;)■• 

M N -2 

. - , £ V 

■tUK 


& = 


(5u( t, x) ) 2 0 (x) dx. 

3 „ 3 , ran 


'n “ “ MD - u (x )) ( A > 

m=ln=l d mn mn 

This energy functional is nonnegative for all 5u and 3<5u/3t in the 

Sobolev space H*(Q) (i.e. the Hilbert space of all real-valued 

square integrable functions defined on Q and vanishing on dQ such that 
their first-order partial derivatives are also square integrable), if 


e V 2 


12EI 


(D - u (x ) ) 3 l 3 A 

d mn mn 


for m = 1, . . . , M, n = 1, . . . , N. 


(28) 


The time rate-of-change of &, after integration by parts, and making 
use of (24) and boundary condition (25), is given by 


M N 


dg/dt = J n ir( I I{( 


c V 2 

o mn 


12EI 


(D - u (x )) 3 l 3 A 

m=ln=l d mn mn 


[Su(t,x ) - <5u(t,x)] 

mn 


e V SV (t) 

o mn mn i , ^ \ i j 

+ }■ 4> (*) dx 

{ \ > 2 I mn 1 

(D - u (x ) ) 

d mn 




(29) 


By requiring that Su(t,x) = Su(t,x ) and 33u(t,x)/3t = 35u(t,x )/3t 

mn mn 

for all x € n , (29) reduces to 


M N 


c V SV (t) (r \ 

de/dt - 1 1 ° “ r ... (l <>*} 

(d - u (x )) v'n * 


m=ln=l 


d mn 


If we set 


35u 


5V (t) = -r 

mn mn Ol 


nn mn 


35u(t,x ) 

mn 

at 


(t,x ) <t> (x), y > 0, 


(30) 


(31) 


then 
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M N 


de/dt - - l l 


e V t 

o mn mn 


(D - u (x ))‘ 
m=ln=l d mn 


/f d> (x ) dxl ( 55u ( t , x )/at) 2 s 0. (32) 

\J n ■" J mn 


The control law (31) implies local rate-feedback at each actuator 
location with feedback gain -y only. 

mn 

Incorporating the foregoing local rate-feedback with the noninter- 
acting controller, the transfer function between u and c is a diagonal 
matrix H(s) whose diagonal elements have the form: 


h n ts) 


g 


(s 2 + 2C us + u 2 ) 
Mi i 


, i = 1 , . . . , NM. 


(33) 


To achieve zero steady-state error for step actuator commands, a 
proportional-plus-integral minor-loop controller is introduced. The 
parameters of the controller and the gain g ( are chosen to ensure 

stability and satisfactory transient response to step-input commands. 
When the noninteracting controls are implemented by means of a digital 
computer, processing delays are introduced. These time-delays may be 
incorporated in h (s). 


5.2 Global Controller for Wavefront Correction 

The initial step in wavefront correction is to estimate the wavefront 
of the incoming light wave reflected from the deformable mirror based on 
the output of the wavefront sensor. In physical situations, a wavefront 
sensor such as the Hartmann-Shack wavefront sensor uses an array of 
micro-lenses, each of which samples a portion of the incoming beam and 
focuses light onto a detector consisting of a CCD camera or a lateral 
field-effect photodiode array. 

Let Q be a bounded open connected subset of R corresponding to the 

s 

effective spatial domain of the wavefront sensor. We introduce a mesh on 

Q whose mesh points are denoted by x mn = ^ X i m ,X 2 n^’ m - 1, . . • n - 

1 N. At each mesh point x , we introduce a patch Q mn , a bounded open 

subset of £1 , representing the effective aperture of the (m,n)-th lenslet 
s 

containing x as an interior point. Let * * »(t,x) denote the wavefront 

mn 

of the incoming wave reflected from the deformable mirror and impinging 
onto the wavefront sensor. The local gradient or angular tilt in the 
wavefront averaged over the aperture of the (m,n)-th lenslet is given by 


V*(t, Ol = f i p (x)V’Ht,x) dx = f _1 (x - x° ), (34) 

x J _mn mn mn 

mn W 

s 

where ^ is a given spatial weighting function associated with the 

mn A q 

(m.n)-th lenslet; f is the focal length of the lenslet; x^ is the 

nominal position of the focal spot for a collimated beam, and x is the 
position of the focal spot for the incoming beam. Thus, the angular tilt 

a o 

in the wavefront can be estimated by measuring the deviations (x - x mn ^- 
From the local gradient data, it is possible to obtain an estimate of the 
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wavefront = 'Ht.x) [12]. For wavefront sensors with circular apertures, 

it is advantageous to express if in terms of Zernike polynomials. By 
comparing if with a reference wavefront 1' r , we obtain the wavefront 

error Sir = if - if. In order to generate the mirror shape correction 
R 

command from 5^, it is necessary to map the wavefront sensor domain I2 g 

onto the mirror domain Q. Let this mapping be a dif feomorphism S from 

onto n. Then the mirror shape correction command 5c is obtained by 

evaluating 5C(t,x) = 5¥(t,S _1 x) at all the actuator locations x^ in ft. 

Again, due to the presence of processing delays in the wavefront 
estimator, a dynamic compensator may be incorporated in the global shape 
controller to ensure overall system stability. 


6. PHYSICAL IMPLEMENTATION 

In the physical implementation of the proposed control system, it is 
desirable to integrate as much as possible the electronic circuitry of 
the controllers with the mirror assembly. Due to the minute capacitances 
associated with the actuators, it is clear that the actuator drivers 
consisting of operational amplifiers must be located in the immediate 
vicinities of the actuators. 

To implement the rate-feedback control given by (31) for damping, 
consider a parallel-plate capacitor driven by an operational amplifier as 
shown in Fig. 5. Assume that the distance A between the capacitor plates 
is time varying so that the capacitance C(t) = c^A/A(t), where A Is the 

area of each plate. Thus, we have 

c (t )ij^2- . = j(t). os) 

where i is the current flowing to the capacitor. Using the expression for 
C(t), (35) can be rewritten as 


c A 

o dv 

"MU "dF 


e o A dA(t) _ 

— V jl = i(t). 


A(t)‘ 


dt 


(36) 


Now, if the voltage v(t) is held at a constant value v^, then the plate 
velocity is related to i(t) by 


dA(t) 

dt 


i(t)A(tr 

c Av 


(37) 


Thus, rate-feedback can be introduced by sensing the current i(t) 
through a resistor in series with the capacitor as shown in Fig. 5. The 
current-sensing resistor can be attached directly to the bottom plate. 

The proportional-plus-integral minor-loop controller can also be 
realized using operational amplifiers which can be integrated with the 
mirror assembly. Since the noninteracting controller requires algebraic 
manipulations, an external digital computer is needed for its 
implementation. However, it is possible to integrate this controller with 
the mirror assembly when single-chip specialized computers become 
available. 

Finally, the wavefront estimator and the global mirror shape 
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controller require a digital computer with sufficiently high speed so 
that the processing time delay will not be detrimental to the performance 
of the overall system. Most likely, these components cannot be integrated 
with the mirror assembly. 

To obtain some information on the orders of magnitude of various 
system parameters and variables, we consider a 1.7 mm square rectangular 
mirror with 15 x 15 actuator patches. Each patch is a 25 pm square pixel. 
The mirror membrane, leaf springs, and supporting posts are micromachined 
from single-crystal silicon sheets. Each leaf spring is a 45-pm long, 
1-pm thick cantilever beam with rectangular cross-section (width = 4 pm). 
Using^ know^data for single-crystal silicon [13^ (Young's modulus E = 1.9 
x 10 N/cm , and mass density p - 2.3 gm/cm ), the electrostatic and 
spring forces associated with a single actuator can be computed from (6) 
and ( 7 ) : 


F = f A = 2.7669 x 10 _21 [V/(D - u)] 2 N 

e e 

(38) 

F = f A = 8.340 x lCf 4 u N 

(39) 


s s 


where V is the actuator voltage, D is in meters, and u is the tip 
displacement of the leaf spring in meters. For a square membrane with 
width L under uniform tension T, the upward restoring force on a square 
pixel with width l located at the center of the membrane with a downward 
displacement u is given approximately by 

F = 8T£u/L N. (40) 

T 

For T = 1 N/m, and the given dimensions for the membrane and pixels, (40) 
becomes 

F = 0. 11765u N. (40' ) 

T 

For a typical membrane displacement u = 0. 1 pm, and D = 0. 5 pm, the 
required actuator voltage V can be computed by balancing the upward 
restoring forces F and F with the downward force F . The resulting V is 

s T e 

equal to 0.8248 volt which is within the operating range of typical 
operational amplifiers. 

Using the foregoing mirror parameters, the static actuator voltages 
V for attaining a bi-parabolic mirror deformation are given by 

mn 

u (x ,x ) = 1.9156 X 10' 13 (850 2 - x 2 ) (850 2 - x 2 ) /im. (41) 

d 1 2 12 

The results are shown in Fig. 6. 

7. CONCLUDING REMARKS 

In this paper , we have only considered the analytical design of a 
control system for a micromachined deformable mirror. The approach is to 
introduce first appropriate static feedback controls for noninteracting 
actuation. Both local rate-feedback and a minor-loop controller are 
introduced for modifying the dyn amics of the actuators. Then a global 
controller is introduced for wavefront correction. Special conside rat ion 
is taken to integrate the controllers with the mirror assembly. Other 
important factors such as thermal effects on the performance of the 
controlled deformable mirror are not studied here. The results pertaining 
to the fabrication of the proposed deformable mirror, and the actual 
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performance of the proposed control system are planned to be reported in 
the near future. 
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APPENDIX 

Explicit expressions for K g and p*^ corresponding to rectangular and 

circular mirrors are given below. ^ 

For a rectangular mirror with spatial domain Q = € ^ : 

i x | < i i = 1 2} and uniform tension T, the Green’s function K is 
i 1 ’ 

given by 


oo oo 4 it ( k7i(x +£ ) ^ ( kn(x +1 U 

Ktx.x.x- x') =1 l TY ~^ — iTT sin ( zH- H 2 1 

s 1’ 2 1 2 L L T(k 2 f 2 + k ,2 £2) l "l 1 \ 2 > 


kn(x +£ )- 
2 2 


wT(kV 

k=l k ’ = 1 1 


( k '* (x ; +£ i^ r * wi *2 t vi 

ln l — 2 T J sin i 2 r J 


k’rr(x’+f )■ 


x s 


(Al) 


Let the actuator patches 0 be square pixels with width A. Thus 0^ - 
{ (x x ) e R 2 ; |x - x | < A/2, |x - x I < A/2). Assuming that the 

1*2 1 lm 2 2n 

actuator weighting function <f> corresponds to the characteristic function 

of U (i e <b (x) = 1 if X e n , and <f> (x) = 0 otherwise), the 

mn ' mn "in ran 

coefficients p lJ defined by (17) are given by 


P IJ = 

mn 


At t <x 8 t , 

1 2 mk nk’ 

z(d - u (x ))\“ / T(k2 < + k,2 0 

d mn k = 1 k -1 1 


oo oo 

1 1 


f k?i(x +£ ) 


x sin 


LX tc K / 

■*»( 


kn(x +£ k 
2J 2l 


2 £ 


where x = (x , x ), an< * 
i j i * J 

fA/2 / k’ n(x’ +f. )> 
sin 


'x -A/2 
lm 


- cos(k’n(x + 

lm 


r +n/2 ( k tux’+£ )> r A « 

* 4 "I* - ) d *; - 

T* *.»}- 


' x -A/2 
2n 


cos(k*7r(x --75“+ i )) 

2n 2 2 


- cos(k’n(x + -«-+£ )) 

2n 2 2 


(A 2 ) 


(A 3 ) 


(A 4 ) 


The Green’s function K for a circular mirror with spatial domain 0. - 

5 

{(r,0), 0 i r s r , 0 s 6 s 2n> and uniform tension T is given by 
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00 00 

} (e, r,0’ ,r’ ) = £ £ 


: ? V ? J (A r) [ J (A r’)cos(k0) 

L 1 irTr 2 [ J (A ,r )] 2 k k “’ tt * 

k=l k ’ =1 o k+1 kk o 


x cos(ke') + J (A r* )sin(k0)sin(k0’ )], (A5) 

kk* kk’ 

where J denotes the Bessel function of the first kind of order k, and 
k 

(A r ) is k’-th zero of J - 

kk' o k 

Let the actuator patches Q be fan-shaped pixels with radial length A 

inn 

and aperture angle 0 as shown in Fig. 3. Thus £2^ = {(0,r): |0 - 0j < 

0 /Z, |r - r I < A/2}. In this case, the coefficients p*^ defined by (17) 

with <p being the characteristic function of £2 are given by 
mn mn 

C 00 00 ry 

5 V V J (A r ) 

mn 2(D - u (x )) 2 A ,, /Tr 2 [J (A r )] 2 k kk J 

d mn k = l k =1 o k + 1 kk o 


A A 

X [a cos (k0 ) + fl sin(k0 ) ] 

ink’ 1 nk’ 1 


where 


+ 

0 /2 


1 ( r r m + 

mk’ 

P 

cos ( k ' 0’ ) 

d0 ’l {I 

i _ 

0 /2 


J \ r - 

mk’ 

P 


m 

1 + 

0 /2 


^ f r r « + 

nk* 

P 

sin(k’ 0’ ) 

de ’} {{ 

1 

0 /2 


J l» r - 

nk’ 

P 


n 


J (A r’ ) r 
kk’ kk’ 


J (X r’ 

kk’ kk’ 


’ dr’ j-. 
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Side view of upper wafer 



Side view of the lower wafer 



Fig. 2 Side view of the deformable mirror. 
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Fig. 5 Implementation of local rate-feedback control. 
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Fig. 6b Required static actuator voltages for bi-parabolic 

mirror deformation given in Fig. 6a. (Maximum voltage at 
mirror center = 0.8248 volt.) 



